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Abstract. Experiments with rare isotopes are shedding light on the role isospin plays 
in the equation of state (EoS) of nuclear matter, and isoscaling -an straight-forward 
comparison of reactions with different isospin- could deliver valuable information about 
it. In this work we test this assertion pragmatically by comparing molecular dynamics 
simulations of isoscaling reactions using different equations of state and looking for 
changes in the isoscaling parameters; to explore the possibility of isoscaling carrying 
information from the hot-and-dense stage of the reaction, we perform our study in 
confined and expanding systems. Our results indicate that indeed isoscaling can help 
us learn about the nuclear EoS, but only in some range of excitation energies. 
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1. Motivation 

Experimental advances in the last decade have allowed the manufacturing of rare 
isotopes -unstable nuclei with excess neutrons- in ways suitable for reactions, this has 
moved the frontier of nuclear science in the dimension of isospin. First results coming 
from studies of rare isotopes shattered the immutability of the shell structure pQ and 
magic numbers and drew a road map to use neutron-rich nuclei to study the elusive drip 
lines which lie far from stable isotopes [2]. These and other results have propelled the 
isotopic degree of freedom to the spotlight, and new facilities and experiments promise to 
help us obtain an improved description of nuclear masses [3], fission barriers, collective 
vibrations [3], neutron skins of neutron-rich nuclei [5] and, in general, understand the 
equation of state (EoS) for neutron-rich matter. 

To reap all these benefits, however, the effect that varying isospin has on the EoS 
must be known and, up to date, this is not the case. Information to constrain the 
nuclear EoS at non-saturation densities and non-zero temperatures has been obtained 
from heavy ion reactions. These collisions produce, albeit for a brief moment, nuclear 
systems at high density and temperature and thus serve as probes of characteristics of the 
sought EoS. This is mainly achieved through comparisons of experimental observables to 
theoretical calculations; this program, however, faces a series of challenges both through 
the detection of clear experimental signatures and the crafting of reliable theoretical 
simulations. 

The main challenge comes from trying to use an EoS to describe a heavy ion 
reaction. In a perfect world, in which the nucleon-nucleon interacting potentials were 
known, excited nuclei in reactions could be described in terms of these interactions. Out 
of these microscopic descriptions of nuclear processes, one could then use appropriate 
averaging procedures to obtain macroscopic information, (e.g. p, e, K, T, p, mean fields, 
enthalpy, etc.) which could be used to craft an EoS. A major condition for this plan to 
work is the attainment of thermal and chemical equilibrium as the macro variables (and 
the EoS) lose meaning in non-equilibrium processes; it is clear that the micro description 
(in a perfect world) would be appropriate to describe all the in- and out-of-equilibrium 
stages of a reaction while the EoS is not. Thus a reaction cannot be fully described 
in terms of the EoS or macro variables, and information about the EoS can only be 
obtained from those segments of the reaction in which the macro variables stabilize; this 
can only be accomplished by means of careful comparison of computational simulations 
of reactions to experimental data. 

To understand the role isospin plays in the EoS, one must first find an observable 
that can provide valuable experimental information that could be analyzed within 
a sensible theoretical framework; initial studies have attempted to do that through 
isoscaling. When comparing the fragmentation yield of stable nuclei reactions with 
those of neutron-rich nuclei collisions, it was observed that they are related by a power 
law of the form Y2(N,Z) oc Yi(N, Z)e~ aN+l3Z , where N and Z are the neutron and 
proton content of the produced fragments, and a and /3 are fitting parameters [6l [7] . 
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The theoretical framework first used to interpret the fitting parameters was that 
of hot and dense fragmenting source disassembling into fragments while equilibrated 
in micro or grand canonical jH], or canonical [9] ensembles. Due to a lack of a 
better approach, the isospin dependence of the binding energy of cold nuclei at normal 
saturation density was assumed to be valid at other densities and higher temperatures to 
find an interpretation of the isoscaling parameters. In this approximation the isoscaling 
parameters were found to be connected to the symmetry energy term of the Weissacker 
mass formula, E sym = C sym (A - 2Z) 2 /A through, a = 4C sym [(Z 1 /A 1 ) 2 - (Z 2 /A 2 ) 2 ]/T, 
where the subscripts refer to the proton and nucleon numbers of the fragmenting sources 
created in the two reactions [10] . 

Unfortunately, the lack of density or temperature dependence of the mass formula, 
as well as other problems, blur the final interpretation of a and (3. A major concern is 
the fact that T and p vary during the reaction, and the concept of a fragmentation source 
is only an approximation useful for invoking statistical models; indeed, the validity of 
C sym (and thus the interpretation of a) have been strongly questioned [IT] . Along the 
same lines, it is known that the isoscaling parameters are know to vary widely during 
the reaction p2]; this casts further doubts about the proposed interpretation of a. 

Thus we arrive at the motivation of the present study: what can we learn about the 
isospin dependence of the EoS from isoscaling? Does isoscaling reflect properties of a 
fragmentation source, or is it set during the expansion of the disassembling source? Can 
isoscaling tell us something about, say, the compressibility of nuclear matter? These 
and other questions will be answered in this study. 

1.1. Plan of action 

Since the process of extracting the isoscaling parameters from experiments is very 
entangled requiring averages and fits over hundreds of different collisions, we feel it would 
be futile to try to connect the parameters to a possible p and T-dependent symmetry 
term, C sym (p,T). Unfortunately, other promising approaches, such as those extending 
the EoS to the non-symmetric case through, e.g., E(p,5) = E(p,5 = 0) + S(p)S 2 , 
where the isospin is inputted through the neutron, proton and nucleon densities, 
8 = (Pn — Pp)J P [E], have not yet been linked to isoscaling. In consequence, we opt to 
follow a much more pragmatic plan of action to understand the connection of isoscaling 
to the EoS. 

In brief, we perform simulations of isoscaling reactions using different equations 
of state and look for any differences in the behavior of the isoscaling parameters; 
same results from different EoS would practically nullify isoscaling as a probe, different 
behavior under different EoS would hold promise. Furthermore, to explore the possibility 
(implied by the proponents of the relationship between a and C sym ) that isoscaling 
carries information from the hot-and-dense stage of the reaction to asymptotia, we 
carry out the proposed study both for systems disassembling under confined conditions 
and during expansion; again, different results would indicate that isoscaling is affected 
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by the expansion, else it would hold promise as a probe of the fragmentation source. 

In section [2] we present argument to justify the computational model used in 
this study, describe how the symmetry energy terms were obtained for each of the 
parameterizations of the potential, and obtain their caloric curves. In section [3] we 
review the concept of symmentropy and study the isoscaling of confined and expanding 
systems. The manuscript closes with a summary of the main conclusions. 

2. Model 

In recent articles we have studied the behavior and properties of isoscaling using both 
classical and geometrical models. In [H] we found isoscaling in simulations of classical 
systems and in [15] we have shown, both analytically and numerically, that isoscaling 
can also be observed in the framework of the nuclear percolation model; an effect totally 
due to the probabilistic aspects of the problem. Along the same line of work, in [16] 
we studied the effect particle correlations have on the isoscaling parameters, and in [17] 
we finally concluded that a minimum isoscaling can be expected for any disassembling 
system based solely on probabilistic aspects. 

In this article we go back to using classical molecular dynamics to study isoscaling 
produced in systems with two different EOS and disassembling in confined and 
expanding environments. In particular, we will use two parameterizations of the 
interaction potential [18] that lead to EoS with compressibility of 250 MeV and 
530 MeV, and will study equilibrated systems clustering in a spherical container as 
well as in a free expansion. 

2.1. Molecular dynamics model 

Here the molecular dynamics model used for this study is introduced along with the two 
potentials that lead to different compressibility values. To study the origin of isoscaling, 
a model capable of reproducing both the out-of-equilibrium and the equilibrium parts of 
a collision is highly desirable; in the present work, we use a molecular dynamics (MD) 
model that can describe non-equilibrium dynamics, hydrodynamic flow and changes 
of phase without adjustable parameters. The combination of this MD code with a 
fragment-recognition algorithm has been applied in recent years to study, among other 
phenomena, neck fragmentation [19] , phase transitions [20] , and other features of nuclear 
reactions, including isoscaling [2T1 [14"] . 

The MD code uses a two-body potential composed of the Coulomb interaction plus 
a nuclear part composed of two Yukawa-like potentials (known as the Illinois potential 
[18J) that correctly reproduces nucleon-nucleon cross sections, as well as the correct 
binding energies and densities of real nuclei. The "nuclear" part of the interaction 
potential is 

V np (r) = V r [exp(-fjL r r)/r - exp(-fi r r c )/r c ] 

- V a [exp(-fi a r)/r - exp(-/j a r a )/r a ] 
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Table 1. Comparison of coefficients obtained for different models 



Coefficient 


Stiff 


Medium 


Experimental 


C v 


16.1 


17.37 


15.75 


C s 


-11.73 


-14.38 


-17.8 


c c 


-0.197 


-0.226 


-0.177 


C sym 


-34.07 


-25.08 


-23.7 



V nn (r) = V pp (r) = Vo [exp(-fi Q r)/r - exp(-fi r c )/r c ] 

where the cutoff radius is r c = 5.4 fm, V np is the potential between a neutron and a 
proton while V nn is that between identical nucleons. The values of the parameters of 
the Yukawa potentials can be selected [T8] as to correspond to an equation of state of 
infinite nuclear matter with an equilibrium density of po — 0.16 fm~ 3 , a binding energy 
E(po) = —16 MeV/nucleon, and a compressibility of ~ 250 MeV for the so called 
medium model, or of ~ 535 MeV for the stiff model. 

2.2. Symmetry energy coefficients 

For completeness, we start by verifying that these two types of interactions indeed 
lead to different values of the symmetry energy for cold nuclei at saturation density. 
To achieve this, we first used dissipative molecular dynamics to construct "nuclei" of 
several masses in their ground states, and then fit the results with the mass formula 
thus obtaining values of the different coefficients, including the symmetry energy term, 

a 

^ sym- 

In a nutshell, in dissipative MD individual nuclei are evolved and cooled until they 
become self-bound and reach proper values of the nuclear radius and binding energy. 
Taking such state as the ground state (still with some residual motion which mimics 
Fermi motion) the operation is repeated for many nuclei of different sizes. We then 
proceeded to fit the resulting nuclei with the liquid-drop mass formula for the nuclear 
binding energy: 

E/A = C V + C S A~V 3 + 4C C (Z {Z - 1)) A~ 4/s + C sym (N - Z) 2 A~ 2 , 

thus obtaining the values of the coefficients presented in Tabled] alongside with published 
experimental data [22] for comparison. In Table[2]we list the nuclei used to fit the binding 
energy, and in Fig. CD we show the binding energies and droplet model values obtained 
for both medium and stiff models. 

These reassuring results indicate that indeed the model used is capable of mimicking 
nuclei with different compressibilities, i.e. with different symmetry energy coefficients. 
Consequently, if isoscaling is to vary under different equations of state, the model most 
likely will show such differences if used with the two parametrizations of the potential. 
We now proceed to use the model to detect differences in isoscaling due to different 
values of the symmetry energy. 
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Table 2. Nuclei used to fit binding energy formula coefficients 



A 


TV 


Z 


EMed 


Estiff 


9 


5 


4 


9.527 


9.617 


12 


6 


6 


10.163 


10.023 


14 


7 


7 


10.148 


10.299 


16 


8 


8 


10.387 


10.376 


20 


10 


10 


10.450 


10.586 


24 


12 


12 


10.521 


10.620 


27 


14 


13 


10.614 


10.801 


35 


18 


17 


10.660 


10.803 


45 


24 


21 


10.580 


10.862 


52 


28 


24 


10.523 


10.819 


59 


32 


27 


10.405 


10.747 


64 


35 


29 


10.360 


10.697 


70 


39 


31 


10.265 


10.618 


79 


45 


34 


10.100 


10.516 


84 


48 


36 


10.026 


10.449 


89 


50 


39 


10.002 


10.395 


96 


54 


42 


9.907 


10.296 


108 


61 


47 


9.744 


10.143 


115 


66 


49 


9.647 


10.054 


119 


69 


50 


9.553 


9.989 



2.3. Simulating the disassembly 

In previous works we have analyzed, among others, the collisions of 40 Ca+ 40 Ca, 
48 Ca+ 48 Ca and 56 Ca+ 56 Ca at different energies, in this work we start by studying the 
evolution of constrained "nuclei" of (JV, Z) = (40,40) and (56,40), which correspond to 
the complete merging of the colliding nuclei of the previous works. 

As done above we have used the dissipative molecular dynamics method to build 
nuclei in their ground states. Once the ground states are available, we place the resulting 
cold systems inside a spherical container. The radii of the container is fixed in order to 
obtain selected values of the number density ranging from 0.001 fm~ 3 to 0.007 fm~ 3 . 
To produce the disassembly, energy is added by scaling the momenta of the particles. 
The trajectories of motion of individual nucleons are then calculated using the standard 
Verlet algorithm with an energy conservation of C(0.01%). We first perform a long 
run in order to let the system relax to equilibrium. Afterwards a much longer run is 
performed and well separated in time snapshots of the evolution are recorded. 

From the microscopic information of the evolution, given by the values of position 
and momenta of the nucleons, we calculate the fragment structure of the system 
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Figure 1. Energies obtained with the mass formula fit (triangles) for the Stiff and 
Medium models (top and bottom panels respectively) together with the corresponding 
ground states calculated using frictional molecular dynamics (circles) 



by means of the MSTE cluster-detection algorithm introduced decades ago [23], but 
recently adapted for the nuclear case [24J. According to this prescription, a particle i 
belongs to a cluster C if there is a particle j in C to which i is bound in the sense 
of < Vij, where pij is the magnitude of the relative momentum, \i the reduced 

mass, and the interparticle potential. In this cluster definition the effect of the 
relative momentum between the particles that form the cluster is taken into account 
in an approximate way. It should be kept in mind that these clusters correspond to a 
given snapshot of the evolution. They are not stable in the sense that their lifetime is 
a function of the particle-particle collision frequency in the system. 

The systems (N,Z) = (40,40) and (56,40), were studied at different energies in 
the range —5 to 8 MeV/A and for four values of the number density with two thousand 
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Figure 2. Caloric curves for the (40,40) system with the medium model for several 
number densities. Squares are for 0.001 /m~ 3 , circles for 0.003 /to -3 , up triangles for 
0.005 /m~ 3 and down triangles for 0.007 /to -3 

snapshots recorded at each energy. It must be remarked that the MD model is fully 
classical and all quantal effects, such as the exclusion principle, Fermi motion and 
isotopic content-modifying phenomena, are excluded. Therefore, any observed variations 
to isoscaling will be entirely due to the change of EOS. [The effect of the Fermi motion, 
although formally absent, is somewhat included by the internal motion of the "nucleons" 
in their "ground state" as explained before.] 

2.3.1. Caloric curves Before proceeding to the study of isoscaling we first characterize 
the disassembling source by its caloric curve, i.e. the functional relationship between the 
temperature and energy of the system. Calculating the temperatures through the usual 
T = (2/3N) Spf/2?n, with p^ representing the momentum of particle i, and the sum 
runs through the ./V particles in the system, the caloric curves for the medium model are 
displayed in Fig. [2] as a function of the energy of the system and for several densitites. 

It is interesting to notice that the systems displays a behavior already found in 
Lennard Jones systems [251 [26], i- e - as the density is lowered the caloric curves start to 
develop a "loop" which signals the presence of a negative specific heat. This behavior 
has been traced to the appearance of "surfaces" in the system, and has been associated 
with a first order phase transition. 

3. Isoscaling 

3.1. Symmentropic limit 

It has been shown that isoscaling is an effect generic to disassemblying systems which 
can be characterized by the concept of symmentropy [T7J [16]; as this is used in the 
following analysis, we first introduce it here. In summary, and for the nuclear case, the 
comparison of two percolating lattices {i = 1,2) with nodes occupied by neutrons and 
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protons with occupation probabilities pz t = Zi/Ai and = Ni/Ai, but similar bond- 
breaking probabilities, produce the ratio of yields R 2 i(N,Z) = Y 2 (N, Z)/Y 1 (N, Z) = 
\pz 2 /Pz^\ Z \Pn 2 /Pn 1 ] N , which is directly related to the isoscaling power law R 2 i(N, Z) oc 
e (aN+/3Z) with a = ln(p N , 2 /p Nl ) and (3 = ln(p Z2 /p Zl ); we take these as the symmentropic 
limit of the isoscaling parameters to be used for comparison. 

As an example, for isoscaling between Ca+ 40 Ca and 48 Ca+ Ca, we have A\ = 80, 
iVi = 40, A 2 = 96, N 2 = 56, and thus p Nl = 40/80, p Ni = 56/96, and a = ln(1.16) = 
0.154; a's near this limit will carry only probabilistic information and will not allow us 
to use isoscaling to probe the EoS. 

3.2. Isoscaling for confined systems 

Data from the evolutions of confined systems as produced in the reactions 40 Ca+ 40 Ca 
and 48 Ca+ 48 Ca were used to construct the yield matrices Yi(N, Z) for the two reactions 
(i = 1, 2) and the corresponding matrix R 2 i(N, Z) = Y 2 (N, Z)/Yi(N, Z). Least squares 
fits to the isoscaling power law yielded values of the parameters a and for each of 
the snapshots of the evolution at the conditions of density and excitation energy stated 
above. 

Fig. [3] shows the behavior of a as a function of excitation energy for different 
densities of the medium model; similar results were obtained for the stiff model. It is 
worth mentioning that a is a decreasing function of the energy of the system, but more 
important, it converges to a constant value with asymptotically appears to approach 
the symmentropic limit; this suggests that at high energies the breakup resembles more 
and more a simple probabilistic disassembly. 

Repeating this type of simulations and analyses for the two interaction potentials, 
can allow us to explore the use of isoscaling to differentiate between a medium- 
compressibility and stiff-compressibility equations of state. For this purpose we have 
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fixed the number density of our confined systems to N/V = 0.007 fm 3 and performed 
the same calculations for both potentials; the resulting values of a are displayed in 
Fig.H 
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Figure 4. a as a function of energy obtained from confined system for the two models 
at a fixed value of the density p = 0.007 fm~ 3 . Squares for the medium model and 
circles for the stiff; the dashed line shows the symmentropic limit. 



Although the medium model has values of a which exceed those of the stiff model 
by as much as ~ 15%, this difference appears to be restricted to the low (self-bound) 
energy regime. At high energies, both compressibilities appear to have values of a that 
are indistinguishable from one another and both converging to the symmentropic limit 
at high energies (see dashed line). 



3.3. Isoscaling for free expanding systems 

To assess the effect of the expansion on isoscaling, we now reproduce the previous 
calculations for systems that underwent an expansion. For this purpose we take the 
configurations recorded for the confined equilibrated systems, remove the confining walls, 
and follow the subsequent evolution by solving the equations of motion. Once that the 
produced fragments and particles are separated and de-excited enough, the steps leading 
to the a parameter are repeated on these asymptotic configurations. 

Fig. [5] shows the values of the a's obtained from these runs. Notice that each point 
in this plot corresponds to a point in Fig. HI and thus allow us to detect a small decrease 
of the overall values of a from confinement to expansion. The asymptotic values of the 
isoscaling parameters indeed are a bit smaller but comparable to those generated at the 
fragmentation stage of the reaction. Therefore, isoscaling appears to have the intrinsic 
ability to provide information of the early moments of the reaction when, presumably, 
the nuclear system is in equilibrium and the information reflects properties of the nuclear 
EoS. 

And that was the good news, the not-so-good- news is that, as found before 
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for confined systems, isoscaling can help us differentiate between stiff and medium 
compressibility EoS only at low energies. Again, systems expanded at high energies 
tend to isoscaling values resembling those predicted by symmentropy and have no EoS 
information to give us. 
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Figure 5. Asymptotic values of a obtained from systems expanding from initial 
densities of p — 0.007 /to~ 3 as a function of their excitation energy. Squares for the 
medium model and circles for the stiff; the dashed line shows the symmentropic limit. 



4. Conclusions 

To quantify the role of isoscaling to provide us with information from the nuclear EoS, 
we have studied isoscaling using molecular dynamic simulations. To assess the relay of 
information from the hot-and-dense stage of the reaction to asymptotia, the study was 
performed with confined systems and compared to similar systems after undergoing an 
expansion. To measure the sensibility of the isocaling parameters to the EoS, the study 
was performed with two systems corresponding to different compressibilities. To further 
connect these calculations with the Weissacker mass formula, the symmetry energy 
coefficient of these medium and stiff media were calculated directly from dissipative 
simulations. 

The disassembling systems studied were equivalent to those formed by collisions 
of 40 Ca+ 40 Ca and 48 Ca+ 48 Ca, and with total binding energies ranging from —4 to 8 
MeV/A. The resulting fragment data were used to construct the ratios R 2 i(N,Z) and 
to obtain the fitting parameters a and (3 of the isoscaling power law; the analysis, 
though, focused on a. 

In summary, studies of both confined and expanded systems and with the different 
interaction potentials demonstrated that: 

i) Confined systems produce slightly larger, but comparable, values of a than 
expanded systems at low energies; the difference disappears at higher energies. 
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ii) In both confined and expanded cases, the medium-compressibility system produces 
slightly larger a than the stiff one, again, this difference is restricted to the low 
excitation energy regime, 

iii) in all cases the values of a decrease with increasing excitation energy converging 
asymptotically to the symmentropic limit in agreement to previous studies [27]. 

These results strongly suggests that: 

iv) isoscaling indeed has the ability to carry information from the early moments of 
the collision, i.e. an asymptotic value of a can be taken as representative of the 
conditions in the hot-and-dense stage of the reaction. 

v) Likewise, the approach to the symmentopic limit, seems to indicate that at high 
energies the breakup approaches a simple probabilistic disassembly; a values at 
these energies do not carry any nuclear information. 

vi) Finally, the use of isoscaling to distinguish between stiff and medium-compressibility 
EoS must be limited to low energies. 

An important suggestion for experimentalists is that, in view of these findings, the 
studies of the symmetry energy through isoscaling must be focused at colliding energies 
leading to relatively cool systems. Another recommendation for future work is the need 
to decouple theoretically the symmentropic and nuclear contributions to the isoscaling 
parameters; we are currently working on such a problem. 

5. Acknowledgments 

C.O.D. is a member of the "Carrera del Investigador" CONICET. C.O.D. acknowledges 
the support of a grant from the Universidad de Buenos Aires, CONICET through grant 
PIP5969. 

References 

[1] P. Rahkila et al, Phys. Rev. C82, 011303 (2010). 

[2] The Frontiers of Nuclear Science, A long Range Plan, The DOE/NSF Nuclear Science Advisory 

Committee, (2007) , www.sc.doe.gov /np/ nsac / docs /Nuclear-Science. Low- Res. pdf . 
[3] P. Danielewicz, R. Lacey, W.G. Lynch, Science 298, 1592 (2002). 
[4] A. Klimkiewicz, et al, Phys. Rev. C76, 051603 (2007). 
[5] B.A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

[6] H. S. Xu, M. B. Tsang, T. X. Liu, X. D. Liu, W. G. Lynch, W. P. Tan, A. Vander Molen, G. 
Verde, A. Wagner, H. F. Xi, C. K. Gelbke, L. Beaulieu, B. Davin, Y. Larochelle, T. Lefort, R. 
T. de Souza, R. Yanez, V. E. Viola, R. J. Charity and L. G. Sobotka, Phys. Rev. Lett. 85, 716 

(2000) . 

[7] M. B. Tsang, W. A. Friedman, C. K. Gelbke, W. G. Lynch, G. Verde and H. Xu, Phys. Rev. Lett. 
86, 5023 (2001). 

[8] M. B. Tsang, C.K. Gelbke, X.D. Liu, W.G. Lynch, W.P. Tan, G. Verde, H.S. Xu, W. A. Friedman, 
R. Donangelo, S. R. Souza, C.B. Das, S. Das Gupta, D. Zhabinsky, Phys. Rev. C64, 054615 

(2001) . 

[9] C. B. Das, S. Das Gupta, X. D. Liu and M. B. Tsang, Phys. Rev., C64, 044608 (2001) 



Dorso, Gimenez Molinelli, Lopez 



13 



[10] A. Ono, P. Danielewicz, W. A. Friedman, W. G. Lynch, and M. B. Tsang, Phys. Rev. C68, 

051601(R) (2003). 
[11] Bao-An Li and Lie- Wen Chen, Phys. Rev. C74, 034610 (2006). 
[12] A. Barran, CO. Dorso and J.A. Lpez, Nuc. Phys. A791, 222 (2007). 
[13] M.B. Tsang, et al., Phys Rev Lett. 102, 122701 (2009). 

[14] C. O. Dorso, C. R. Escudero, M. Ison and J. A. Lopez, Phys. Rev. C 73, 044601 (2006). 
[15] A. Davila, C. Escudero, J. A. Lopez and C. O. Dorso, Physica A 374, 663 (2007). 
[16] C. O. Dorso, Phys. Rev. C73 034605 (2006). 

[17] L. Moretto, C. O. Dorso, L. Phair and J.B.Elliott, Phys. Rev. C77 037603 (2008). 

[18] R. J. Lenk, T. J. Schlagcl and V. R. Pandharipande, Phys. Rev. C42, 372 (1990) 

[19] A. Chernomoretz, L. Gingras, Y. Larochelle, L. Beaulieu, R. Roy, C. St-Pierre and C. O. Dorso, 

Phys. Rev. C65, 054613 (2002). 
[20] A. Barrafion, C. O. Dorso and J. A. Lopez, Rev. Mex. Ffs. 47-Sup.2, 93 (2001). 
[21] A. Barrafion, C. O. Dorso, and J. A. Lopez, Nuclear Phys. A791, 222 (2007). 
[22] J.M. Rohlf, Modern physics from A to Z, John Wiley and Sons, (1994). 
[23] T.L. Hill, J. Chem. Phys 23, 617 (1955). 
[24] A. Strachan and C. O. Dorso, Phys. Rev. C56, 995 (1997). 

[25] A.Chcrnomoretz, M. Ison, S. Ortiz and C.O.Dorso , Phys. Rev. C64, 024606 (2001). 
[26] M. J. Ison and C.O.Dorso, Phys. Rev. C69, 027001 (2004). 

[27] C. O. Dorso, C. M. Hernandez, J. A. Lopez and J. A. Mufioz, Phys. Rev. C78, 034613 (2008) 



